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We present results on the static three- and four-quark potentials in SU(3) and 
SU(4) respectively within quenched lattice QCD. We use an analytic multi-hit pro- 
cedure for the time links and a variational approach to determine the ground state. 
The three- and four-quark potentials extracted are consistent with a sum of two-body 
potentials, possibly with a weak many-body component. The results give support 
to the A ansatz for the baryonic area law. 
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The nature of the three quark potential is of prime importance in the understanding of baryon 
I—I . structure. However up to now it has received little attention in lattice QCD studies. This is to be 
^v ' contrasted with the quark - antiquark potential relevant for meson structure for which many lattice 
results exist Q. 

The aim of the present work is to study the nature of the three quark potential within lattice 
QCD. The fundamental question which has been raised more than twenty years ago, is whether the 
static three quark potential can be approximated by a sum of three two body potentials, known 
in the literature as the A-ansatz, or whether it is a genuine three body potential. The latter is 
obtained in the strong coupling approximation by minimization of the energy of the three quark 
state. The resulting minimal length flux tube is a configuration where the flux tubes from each 
f~^ ■ quark merge at a point. Due to its shape it is known as the Y"-ansatz. 

Recently two lattice studies of the three quark potential have reached different conclusions: Pre- 
liminary results by G. Bali B at (3 = 6.0 favour the A ansatz whereas the analysis of lattice results 
at (3 = 5.7 by Takahashi et al. |2| gives more support to the Y— ansatz. The difficulty to resolve 
the dominant area law for the baryonic potential is due to the fact that the maximal difference 
r"| , between the two ansatze is a mere 15%. 

In our study we make a number of technical improvements in order to try and distinguish the Y — 
and A— ansatze. Besides using the standard techniques of smearing and the multi-hit procedure for 
noise reduction, we employ a variational approach y to extract the ground and first excited state 
of the three quarks. Both the multi-hit procedure, which is done analytically, and the variational 
approach were not used in ref. B] . These are especially important for the larger Wilson loops where 
the confining part of the potential is the most dominant. Instead of the multi-hit procedure for the 
time links we have also tried the recently proposed hypercubic blocking H. We did not, however, 
find any improvement as compared to the multi-hit procedure. 

In addition to the SU(3) gauge group we also present results for SU(4). Since the same issue 
of which ansatz is favoured arises in any gauge group a calculation in SU(4) can help decide the 
preferred area law. The difference between the two-body approximation and the many-body force 
is bigger for SU(4), reaching for the lattice geometries that we looked at a maximum value of 20%. 

The SU(3) baryon Wilson loop is constructed by creating a gauge invariant three quark state at 
time t = which is annihilated at a later time T. 

W 3q = ±e abc e a ' h ' c 'U{x,y,l) aa 'U(x,y,2) bb 'U(x,y,3r' (1.1) 
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for the three quark lines that are created at x and annihilated at y and 

U(x,y,j) =Pexp ig I dx^A^x) 



(1.2) 



where P is the path ordering and T(j) denotes the path from x to y for quark line j as shown in 
Fig. |. 
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FIG. 1. The baryonic Wilson loop. The quarks are located at positions ri,r*2 and 7-3. 

The three quark potential is then extracted in the standard way from the long time behaviour of 
the Wilson loop: 



V 3q = - lim - In < W 3q > 



(1.3) 



In SU(4) the corresponding colour singlet gauge invariant four quark state is constructed in an 
analogous manner and the four quark potential Vi q is similarly extracted. We will be using the 
term baryonic potential to denote the colour singlet combination of N quarks despite the fact that 
in SU(4) the spin of the four quark state is an integer. 



2. WILSON LOOP FOR THREE AND FOUR QUARKS 



2.1. SU(3) 



We describe in more detail here the two possibilities put forward for the area-law of the SU(3) 
baryon Wilson loop. In the strong coupling limit in the presence of three heavy quarks the gauge 
invariant three quark state with the least amount of flux will yield the lowest energy. If the three 
quarks are at positions ri , V2 and r3 and provided none of the interior angles of the triangle with 
vertices at the quark sites is greater than 120° then the flux tubes from the quarks will meet at an 
interior point r$ || . The position r^ is determined by minimizing the static energy with the result 
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which is known as the Steiner point. The angles between the flux tubes are 120° independently of 
the vectors r^. If one of the interior angles of the triangle of the quarks is greater than 120° then the 
flux tube at that angle collapses to a point. Time evolution of this state produces a three-bladed 
area. This area law is the Y-ansatz mentioned in the Introduction. We denote the minimal length 
of the flux tube for this ansatz Ly and the corresponding area Ay . 

The second possibility is that the relevant area dependence of the baryonic Wilson loop is 
given by the sum of the minimal areas A^ spanning quark lines i and j. This is known as the 
A-ansatz with the corresponding length and area denoted by L& and A& respectively. 

The position of the Steiner point can be obtained analytically H in terms of the three quark 
positions and the difference between the two laws as compared to the two-body ansatz, 

J2 r i*- n^2 r Jk)/^^2r jk , (2.5) 



2 , _ 

j<k ' j<k 

attains @J the maximum value of (Ly — £a/2)/(£a/2) = 2/sqrt3 — 1 = 0.154.. when the quarks 
form an equilateral triangle. The factor of 1/2 is due to the non-Abelian nature of the gauge 
couplings giving half as much attraction for a qq in an antisymmetric colour state as a qq in a 
colour singlet. In general the attraction for (N — 1) quarks in an N quark antisymmetric colour 
state is a factor 1/(JV— 1) less that the attraction for a qq in a colour singlet. Because of this factor 
W(iVi) < Ly. 



2.2. SU(4) 

In SU(4) the ground state of the system in strong coupling corresponds to the configuration with 
minimal length for the flux tubes which join the quarks. Minimization of the static energy results in 
the introduction of two Steiner points, A and B somewhere in space, with the flux tubes from two 
quarks joining at A, while the flux tubes from the other two quarks meet at B. This configuration 
is visualized in Fig. @. Since 4x4 = 6010 the two lines emanating from the two Steiner points 
join to form a colour singlet. In analogy to SU(3) we will call this area law as the y-ansatz, with 
a corresponding flux tube length Ly. 

The two Steiner configuration is always favoured as compared to a single Steiner point defined 
by equation 

Y p ~ rA \ = o . (2.6) 

r-f \n -r A \ 
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where all four quark lines meet and which is a local minimum of the static energy. The area law for 
the baryonic Wilson loop takes now the shape of a long four-bladed surface with the blades meeting 
at A as shown in fig. 0. Due to this shape, we refer to this configuration as the A-law and denote 
the corresponding flux tube length as Lx- 

In contrast to SU(3) where for any given location of the three quarks, the Steiner point and 
therefore the energy can be computed analytically, in SU(4), the two Steiner points in the F-ansatz 
can be obtained by a simple iterative numerical procedure. The two Steiner points have vectors 
that meet each at 120° and one Steiner point can be obtained in terms of the other. Starting from 
an initial guess for the position of one of the Steiner points, r^, we can compute rs as the Steiner 
point of r 3 , r 4 , ta- The r%,r2, rg vectors lead now to a new estimate for a Steiner point r A which 
in turn is used to compute a new Tb etc. [The procedure converges after 30-40 iterations to the 
minimum.] The location of the single Steiner point is easily computed by a numerical solution to 
eq. 2.6. 






FIG. 2. The Wilson loop for four quarks. The quarks are located at positions 7-1,7-2, 7-3 and 7*4. The 
upper graph shows the local minimum of the energy with one Steiner point A, and the lower is the minimum 
with two Steiner points A and B. 

It has been argued in [J7| that the two-body force is the relevant interaction for any SU(N) 
gauge theory. It is proven in |7| that Ly > La/(N — 1) holds for any location of the four 
quarks. From the numerical investigation , it turns out that the relative difference between the Y— 
energy and the two-body law is maximal for the configuration of maximal symmetry for the four 
quarks. This amounts to putting the quarks on the vertices of the regular tetrahedron and gives a 
relative difference of 21.96 % with respect to the two-body term. This difference decreases as the 
configuration becomes more asymmetric in space and can decrease down to 5 — 6 % for the most 
asymmetric locations of the quarks on a 16 3 lattice. Therefore, in order to obtain a clear signal on 
which law is preferred by the £{7(4) quarks, we studied geometries with maximal symmetry. 

As far as the four-bladed surface area law is concerned, we observed that Lx always exceeds Ly 
by at most 3.5 %. In fact, the ratio, (Lx — Ly)/ Ly, becomes minimal for the most symmetric 
configuration of the tetrahedron, obtaining a value of just 0.43 %. Here the ratio in fact increases 
as one increases the asymmetry of the four quark locations, becoming maximal if all four quarks are 
located on a plane. In particular, if the quarks are located at the vertices of a square, (Lx — Ly)/Ly 
takes its maximal value of 3.5 %. With the current data, discriminating an effect of 0(3) % between 
the Y- ansatz and the X— ansatz is not possible. Therefore, we will pick geometries that maximize 
the difference between the Y- and A ansatze. Since for these geometries the difference between the 
Y- and X-ansatze is of the order of 2% one has to keep in mind that when we refer to the F-ansatz 
we will in fact mean the area law with one or two Steiner points. 

The factor of 1/(N — 1) which relates the long range part of the two-body qq and qq potentials 
also occurs in lowest order gluon exchange so that the two-body short range potential is given by M 
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where Vkj is the qq one-gluon potential 



V jk = - 
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with C F = (N 2 - 1)/2N the quark Casimir of 0(N). 

Thus the expected forms of the "baryonic" potential in SU(N) that we will be applying to SU(3) 
and SU(4) are 

W r 1 ,., r „) = ^ - F l TS |^ + ¥ l T . i4 (2.9) 
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with a the string tension of the qq potential. 



, -rjk 



3. LATTICE TECHNIQUES 

As we have mentioned in the Introduction, the two recent lattice studies of the baryonic potential 
|g,y have yielded different conclusions, the first supporting the Y— ansatz and the second the A— 
ansatz. Since the difference between the two ansatze is ~ 15% for SU(3), obtaining conclusive 
results requires making a large effort to reduce the statistical noise, especially for the large loops 
where the absolute difference between the two ansatze becomes more visible. In this work, we used 
a number of improvements as compared to previous studies in SU(3). To our knowledge, this is the 
first measurement of the 4-quark potential in SU(4). We describe briefly the techniques that we 
use in order to reduce noise and extract more reliably the ground state. 
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We use the multi-hit procedure g] for the time links. For SU(3) the temporal links are 
integrated out analytically @] and substituted by their average value 

f AT] TJ.( n ) e PSi{U) 

with Si(U) = jjTr (1/4(11,) F^ (n)) and F(n) is the staple attached to the time link that is being 
integrated over. It has been shown in SU(2) [|| that replacing the time links by their average 
value in this fashion reduces the error on large Wilson loops of the order of tenfold. The 
factor found in ref. is x 2T ~ 0.889 2T where T is the time extent of the Wilson loop. For 
the SU(3) baryon loop the reduction factor will be x 3T giving an even larger noise reduction 
for the large loops. We point out here that the multi-hit procedure was not used in ref. M. 
In SU(4) the integration over the temporal links was done numerically. 

We compared the multi-hit procedure with the recently proposed hypercubic blocking M 
on the time links. Using the optimal parameters given in [H at /3 = 6.0, we compare in 
Fig. ra the results on the same configurations, using the analytic multi-hit procedure and 
using hypercubic blocking. As it can be seen, the multi-hit procedure gives smaller errors for 
large loops and therefore we adopt it in this work. 
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FIG. 3. The qq potential for SU(3) at /3 — 6.0 on a 16 3 x 32 lattice using the multi-hit procedure with 
x-symbol and hypercubic blocking shown with the crosses. 



• To maximize the overlap of the trial state with the three quark ground state we use smearing 
of the spatial links || . We replace each spatial link by a fat link by acting on it with the 
smearing operator S defined by 



SUj(x) =v(Uj(x) + aJ2 [Uk(x)Uj(x + ak)u£(x + aj) 



(3.12) 



where V denotes the projection onto SU(3). This is iterated n times. We consider M different 
levels of smearing and construct anMxM correlation matrix of Wilson loops |1C]] . For the 
parameter a and the number of smearings, ni, for each different smearing level I we take what 
is found to be optimal in p0[; namely 
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for smearing levels / — 0, • • • , M — 1 and r$ Sommer's reference scale [jllj . In all our computa- 
tions we used M = 4. For SU(4) at j3 = 10.9 we found that the parameters used for (3 = 5.8 
in SU(3) produce reasonable results. 

The correlation matrices C{t) for the mesonic and baryonic Wilson loops were analyzed using 
a variational method H. We use two different variants both yielding consistent results. 



In both variants we solve the generalized eigenvalue problem |10 

C{t)v k {t) = X k (t)C(to)v k (t) 

taking to/a = 1. In the first variant, the potential levels are extracted via 
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by fitting to the plateau. In the second variant we consider the projected Wilson loops 



W P (t) = v^(t )C(t)v (to) 



(3.16) 



and fit to the plateau value of — In I Wp(t+ \)/Wp{t) I . In Fig. [| we show the results of these 
two variants for SU(4) for the four quark static potential. 
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FIG. 4. Comparison of the plateaus obtained in SU(4) for — ln(Ao(£ + l)/Ao(t) solving the generalized 
eigenvalue equation at each time t (pluses) and with the projected Wilson loop —In I Wp(t + 1)/Wp(t) 
(x-symbols). 

The projected correlation has a larger contamination of excited states for time slice t/a = 3 
but by the next time slice the two procedures yield the same results. We have found that 
for SU(3), where the ni for the l-th smearing level are larger for the corresponding j3 values 
than the number of smearings used in SU(4), the projected method yields smaller errors. In 
all cases we have checked that the values we extract for the ground state within these two 
procedure are consistent with each other. From Eq.( (3.15 ) we also obtain the energy for the 
first excited state. Although the data are rather noisy, we can obtain an estimate, which we 
use to fix the minimum value for the time interval used in the extraction of the ground state, 
such that the excited state contamination is less than e -2 . 



4. RESULTS 

All the computations were carried out on lattices of size 16 3 x 32 at (3 values 5.8 and 6.0 for 
SU(3) and 10.9 for SU(4). The (3 value for SU(4) was chosen so that the lattice spacing is close 
to the value for SU(3) at (3 = 6.0. In the case of SU(3) we used 200 configurations at (3 = 5.8 
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and 220 at j3 = 6.0 available at the NERSC archive pjj and for SU(4) we generated 100 quenched 
configurations. 

We consider geometries on the lattice which produce the biggest difference between the A— and 
Y *— ansatze. For SU(3) each quark is placed on a different spatial axis equidistant from the origin. 
The results are shown in Figs. H and y for /3 — 5.8 and 6.0 respectively. To reduce systematic 
errors when comparing with the qq potential, we also compute, on the same configurations, the 
static qq potential with the quark and the antiquark at the same locations as the 3 quarks of the 
qqq potential. The errors shown on these figures are the jackknife errors. The string tension in 
lattice units extracted from fitting the qq potentials is a^fa = 0.329(3) at (3 = 5.8 and 0.224(3) 
at (3 = 6.0 consistent with the value of ref. p3|| . At short distances the baryonic potential, V 3q , is 
approximately equal to the sum of the corresponding two - body potentials i.e. we find agreement 
with the tree level result that V 3q ~ 3/2V 9 q. At larger distances, V 3q is enhanced compared to the 
tree level result. On the same figures we also show the curves corresponding to the A— and Y— 
ansatze. The lattice data lie closer to the curve given by the A— area law. However, at distances 
larger than about 0.7 fm, the three-quark potential appears enhanced as compared to the sum 
of the two-body potentials. This enhancement can be explained by a small admixture of a three 
body-force, although it is is so small that it might also reflect imperfections in our variational search 
for the groundstate. 

In SU(4) we studied three different geometries chosen so that the difference between the A and 
analog of the Y law is maximal. In what we call geometry 1 the quarks are placed symmetrically 
on a plane distance / from the origin. The energy difference between the two ansatze is 20.0%. In 
geometries 2 and 3, three quarks have coordinates (1,0,0), (0,1,0), (0,0,1) whereas the fourth is at 
(0,0,-1) for geometry 2 and at the origin for geometry 3. The energy differences between the A and 
Y laws are 20.1% for geometry 2 and 19,1% for geometry 3. The string tension is obtained by 
fitting the on axis qq potential excluding the first point. We find aa = 0.238(4) in agreement with 
the value of 0.2429(14) of ref. @]. The quality of the fit is shown in Fig. | with x 2 /d-cf = 1.0 
where we included the results when the quark and the antiquark are on different axes. 
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FIG. 5. The static baryonic potential at (3 — 5.8 (pluses) in physical units. The crosses is the sum of 

the static qq potential. The curves for the A and Y ansatze are also displayed. The quarks are located at 
{1, 0, 0), (0, 1, 0), (0, 0, 1) and r = r 12 = n 3 = r 2S = V^l . 
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FIG. 6. As figure | but for (3 - 6.0 
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FIG. 7. The <?g potential for SU(4) at f3 — 10.9 fitted to the form Vo — b/r + ar. The errors shown are 

the jackknife errors. 

The corresponding results for the four-quark static potential are shown in Figs, g, || and [h] for 
geometries 1, 2 and 3 respectively. Again we find that the four quark potential is approximated by 
the sum of qq potentials with a small enhancement at larger distances. The results in all cases lie 
closer to the A ansatz. 
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FIG. 8. As figure g but for SU(4) geometry 1. The quarks are located at (Z,0,0), (0,1,0), (-1,0,0), 

(0, —I, 0) and r = n 2 = r 2 3 = r 34 = ru = V% I- 
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FIG. 9. As figure | but for SU(4) geometry 2. Here the quarks are located at (1,0,0), (0,1,0), (-1,0,0), 
(0, 0, 1) and r = r 12 - r 2 3 = r 34 , = r i4 = y/2 I. 
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FIG. 10. As figure I but for SU(4) geometry 3. The quarks are located at (r,0,0), (0,r, 0), (0,0, r) and 
(0,0,0). 
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5. CONCLUSIONS 

Our results for the static three- and four-quark potential in SU(3) and SU(4) are consistent with 
the sum of two-body potentials below a distance of about 0.8 fm, and clearly inconsistent with the 
Y— ansatz. 

For larger distances, where our statistical and systematic errors both become appreciable, there 
appears to be a small enhancement due to an admixture of a many-body component. Nevertheless, 
for the distances up to 1.2 fm that we were able to probe in this work, the A area law gives the 
closest description of our data. 

We have made use of all the known techniques in order to reliably identify the plateaus in the 
Wilson loops and extract the ground state potential. Nevertheless, for the larger loops the plateaus 
were hard to identify, resulting in large errors. This is a challenging numerical problem, and we 
cannot exclude the possibility that the small enhancement of the potential above the A area law 
which we observe is simply caused by a failure to filter out all excited states in our variational 
search for the groundstate. Taking the results in both SU(3) and SU(4) at face value the conclusion 
that can be drawn is that the A area law provides the closest description to the baryonic potential 
up to distances of 1.2 fm. More refined techniques for noise reduction for the large loops will be 
needed in order to clarify whether a genuine many-body component is present at larger distances. 
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